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ABSTRACT 

Large synoptic (repeated scan) imaging sky surveys are poised to observe enormous 
numbers of core-collapse supernovae. We quantify the discovery potential of such sur- 
veys, and apply our results to upcoming projects, including DES, Pan-STARRS, and 
LSST. The latter two will harvest core-collapse supernovae in numbers orders of mag- 
nitude greater than have ever been observed to date. These surveys will map out the 
cosmic core-collapse supernova redshift distribution via direct counting, with very small 
statistical uncertainties out to a redshift depth which is a strong function of the survey 
limiting magnitude. This supernova redshift history encodes rich information about 
cosmology, star formation, and supernova astrophysics and phenomenology; the large 
statistics of the supernova sample will be crucial to disentangle possible degeneracies 
among these issues. For example, the cosmic supernova rate can be measured to high 
precision out to z ~ 0.5 for all core-collapse types, and out to redshift z ^ 1 for Type 
Iln events if their intrinsic properties remain the same as those measured locally. A 
precision knowledge of the cosmic supernova rate would remove the cosmological un- 
certainties in the study of the wealth of observable properties of the cosmic supernova 
populations and their evolution with environment and redshift. Because of the tight link 
between supernovae and star formation, synoptic sky surveys will also provide precision 
measurements of the normalization and z < 1 history of cosmic star-formation rate in 
a manner independent of and complementary to than current data based on UV and 
other proxies for massive star formation. Furthermore, Type II supernovae can serve as 
distance indicators and would independently cross-check Type la distances measured in 
the same surveys. Arguably the largest and least-controlled uncertainty in all of these 
efforts comes from the poorly-understood evolution of dust obscuration of supernovae 
in their host galaxies; we outline a strategy to determine empirically the obscuration 
properties by leveraging the large supernova samples over a broad range of redshift. 
We conclude with recommendations on how best to use (and to tailor) these galaxy 
surveys to fully extract unique new probes on the physics, astrophysics, and cosmology 
of core-collapse explosions. 
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Introduction 



A new generation of deep, large-area, synoptic (repeated-scan) galaxy surveys is coming online 
and is poised to revolutionize cosmology in particular and astrophysics in general. The scanning 
nature of these surveys will open the way for a systematic study of the celestial sphere in the 
time domain. In particular, ongoing and planned surveys are sensitive to the transient cosmos on 
timescales from hours to years, and to supernova flux limits down to 24^^^^ and sometimes fainter. 
As we will see, these capabilities will reap a huge harvest in cosmic supernovae and will offer a new 
and direct probe of the cosmic supernova history out to high redshifts. 

In the past decade, supernovae in nearby and distant galaxies have come to pl ay crucial role 



for cosmology, v ia the use of Type la explosions as "standardizable" candles (e.g., iPhillips 



1993 



Riess et al.lll996l ). These powerful beacons are detectable out to very high redshift and thus reveal 
the cosmic expansion history for much of the lifetime of the universe; the stunning result has been 
the detection of the acceleration of the Universe and the infer ence that dark en e rgy of some form 



dominates the mass-ene r gy content of the cosmos today (e.g., iRiess et al.lll998l : iPerlmutter et al 
19991 : lAstier et al.l l2006l : IWood-Vasey et al.l 120071 ). The detection of large numbers of Type la 



supernovae over a large redshift ra nge, and their use as cosmological probes, represents a major 
focus of future galaxy surveys (e.g.. lWang et al.ll2004l ). 



While studies of supernova Type la (thermonuclear explosions) justly receive enormous atten- 
tion due to their cosmological importance, there has been relatively little focus on the detection of 
the more numerous population of core-collapse supernovae. These explosions of massive stars show 
great diversity in their observed properties, e.g. including several varieties of Type II events, but also 
Types lb and Ic events. Despite their heterogeneous nature, some core-collapse events may nonethe- 
less provide standardized cand les, via their early lightcurves whose nature is set by the physics of 
their expanding photospheres (jKirshner &: KwanI Il974l : I Baron et al.l |200J; iDessart &: Hillieii |2005| . 
see below). Moreover, core-collapse events are of great intrinsic importance for cosmology, astro- 
physics, and particle physics. These events play a crucial role in cosmic energy feedback processes 
and thus in the formation and evolution of galaxies and of cosmological structure. 

Synoptic surveys tuned for Type la events will also automatically detect core-collapse super- 
novae. Indeed, as survey coverage and depth increase, they will, for the first time, image a large 
fraction of all unobscured cosmic core-collapse supernovae out to moderate redshift. These pho- 
tometric detections of supernovae and their light curves will shed new light on a wide variety of 
problems spanning cosmology, particle astrophysics, and supernova studies. Moreover, these data 
will "come for free" so long as surveys include core-collapse events in their analysis pipelines. 



For example. iMadau et al.l (|l998l ) already pointed out the link between the cosmic star forma- 
tion history and the cosmic supernova history, and showed that when integrated over all redshifts, 
the all-sky supernova event rate is enormous, ~ 5 — 15 events/sec in their estimate. Upcoming 
synoptic surveys will probe most or all of the sky at great depth, and thus are positioned to ob- 
serve a large fraction of these events. Consequently, these surveys will reveal the history of cosmic 
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supernovae via directly counting their numbers as a function of redshift. 



Already, recent and ongoing surveys have begun to detect core-collapse supernovae. However, 
to date, surveys have focused on Type la events, and thus core-collapse discovery and observation 
has been a serendipitous or even accidental byproduct of SNIa searches. As a result of these surveys, 
the supernova discovery rate is accelerating, and the current all-time, all- Type supernova count is 
~ 5000 since SN1006l^ Thus core-collapse data is currently sparsely analyzed and reported in an 
uneven manner. This situation will drastically improve in the near future, when the supernova count 
increase by large factors, culminating in up to ~ 100,000 core-collapse events seen by LSST annually. 
In this paper we therefore will anticipate this future, rather than make extensive comparison with 
the present data though we will make quantitative contact with current results. 



Our work draws upon several key analyses. The thorough and elegant work of lDahlen &: Fransson 



(jl999l ) laid out the framework for rates and observability of cosmic supernovae of all types. Their 
work assembled a large body of supernova data and applied it to make rate and discovery predic- 
tions for the wide variety of star formation histories and normalizations viable at that time, with 
a particular focus on forecast s for very high redsh ift (out to z ~ 5) observable by the infrared 
James Webb Space Telescope. I Sullivan et al.l (|200d ) estimated the rates for supernovae lensed by 
the matter distribution-particular ly rich clusters-along the line of sight; these objects further ex- 
tend the reach of infrared searches, and identified a possible su pernova candida t e from Hubble Space 
Telescope archival images of an intermediate-redshift cluster. I Gal- Yam et al.l (120021 ) made similar 
calcul ations of the infrar e d obs ervabi lity of supernovae, an d identified additional events in archival 
data. iGal-Yam Sz Maoa (120041 ) and lOda &: Totanil (j2005l ) presented forecasts for then-upcoming 
ground-based surveys. These studies considered all cosmic supernovae, but with a focus on Type la 
events, specifically with an eye towards revealing the Type la delay time as well as a parameterized 
characterization of the cosmic star formation history based on Type la counts. In addition, these 
first studies reasonably chose to emphasize near-term (i.e., now-completed or ongoing) relatively 
modest surveys, or on future space-based missions such as SNAP, with little to no study of the 
impact of large synoptic surveys. Moreover, while these works included dust extinction effects in 
host galaxies, but because of th eir focus on Type la e vents, they did not study the possibility of a 
redshift evolution in extinction (jMannucci et al.l 120071 ). 



We buil d on the important st udies of lDahlen Sz FranssonI (j 19991 ) , iGal-Yam et al.l (|2002l ) ; iGal-Yam Sz Maoz 



(|2004l ). and lOda &: Totanil (j2005l ) in several respects: (1) we explore the promise of synoptic sur- 
veys and forecast the very large numbers of supernovae they will find; (2) we focus on less-studied 
core-collapse ev e nts; (3) we incorporate the (pessimistic) possibility of strong dust evolution of 
Mannucci et al.l (120071 ) which is a dominant obstacle to observing massive star death at high red- 
shifts; (4) we present a strategy for empirically calibrating the obscuration properties across a broad 
range of redshift by studying the evolution of the supernova luminosity function; and (5) we study 
the unique opportunities that become available with the large supernova harvest of synoptic sur- 



^ Central Bureau for Astronomical Telegramsl 1 20081 ): see also |http : //www . cf a . harvard . edu/iau/lists/Supernovae . html] 
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veys; in particular, we show how the cosmic supernova rate can be recovered based on core-cohapse 
counts, without assumption as to its functional form. 

Our goal in this paper is to explore the impact synopic surveys will make on core-collapse 
supernova astrophysics and cosmology. We summarize key upcoming surveys in ^ In ^ we 
review expectations for the CSNR, core-collapse supernova observables, and the effect of cosmic 
dust and its evolution. We combine these inputs in [2] where we forecast the core-collapse supernova 
discovery potential for upcoming surveys. We quantify in detail the strong dependence of the 
supernova harvest on the survey limiting magnitude, which we find to be the key figure of merit 
for supernova studies. We discuss some of the supernova science payoff in ^ and conclude in ^ 
with some recommendations for synoptic surveys. 



2. Synoptic Surveys 



Current ar id future sky su rveys build on the pioneering approach of the Sloan Digital Sky 
Survey (SPSS: lYork et al.ll200Cl ). Following SDSS, these surveys will produce high-quality digital 
photometric maps of large regions of the celestial sphere. The powerful innovation the new surveys 
to extend the original SDSS approach into the time domain. Each program will scan part of their 
survey domain frequently, with revisit periods of days and in some cases even hours, and maintain 
this systematic effort throughout the survey's multi-year operating lifespans. The result will be 
unprecedented catalogs of transient phenomena over timescales from hours to years. These surveys 
are thus ideal for supernova discovery and matched to supernova light curve evolution timescales; 
the result will be a revolution in our observational understanding of supernovae. 



Table 1: Recent and Future Synoptic Sky Surveys 



Survey 


Scan Area 


SN Depth 


Scan 


Expected 


Name 


AJlscan [deg2] 


r-band mf^ [mag] 


Region 


Operation 


SDSS-II 


300 


21.5 


SDSS southern equatorial strip 


2005-2008 


DES 


40 


24.2 


South Galactic Cap 


2011-2016 


Pan-STARRS 


30000 


23 


~ 75% of the Hawai'ian sky 


2010-2020 


LSST 


20000 


23-25 


southern hemisphere 


2014-2024 



The science harvest in the time domain depends on both the depth of the scans and their 
breadth across the celestial sphere. These scale with collecting area A and sky coverage l^survey, 
respectively. Consequently, the figure of merit for scan power is the etendue Ar^survey Forthcoming 
projects are designed to maximize this quantity. 

The viability of supernova discovery, typing, a nd followup by large- scale synoptic surveys has 
now b een tested by the SDSS-II supernova search (jFrieman et al.l 120081 ). This program extended 
SDSS ( York et ahlboOC ) into the time domain, scanning at a ~ 5 day cadence, identifying and typing 
supernova candidates from photometric data in real time, and following up with spectroscopic 
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confirmation. This survey will serve as a testbed for the larger future campaigns. It is thus very 
important and encouraging that SDS S-II has reported the discovery of 403 confirmed supernovae 



in the first two seasons of operation (jSako et al 



20081 ). The search algorithms and followup were 



focused on Type la events, for which light curves and spectra have been recovered over 0.05 < 
z < 0.35; human input was used for supernova typing, but aut omated routines a.p pear promising 
and will be essential for larger surveys. Follow-up spectroscopy (jZheng et al.ll2008l ) yields accurate 
supernova and host-galaxy redshifts (u™ ^ 0.005 and crP ^ 0.0005); host-galaxy contamination is 
found to be well-addressed by fitting and a principal component analysis. 



Table [T] lists several major current and future synoptic surveys, and gives the values or current 
estimates of their performance characteristics. The mf^ values are derived from the survey 5a 
detection s for single visit exp osures, which have been corrected l™'^^ shallower as noted above. 



SDSS-II (jFrieman et al.ll2008l ) is recently completed, as discussed above; we adopt an r-band lim- 



iting magnitude of 21.5™'^^ (J. Frieman, priva te communication). The Dark Energy Survey (DES; 
The Dark Energy Survey Collaboration! 120051 ) will push down to nii^ ~ 24. 2™'^^ in r-band; as 
we will see below, this will already enormously incr ease the sup e rnova harves t. Finally, look- 
i ng out farther into the next decade, Pan- STARRS (| JewittI l2003l : iTonrvl |2003| ) and then LSST 



( The LSST Collaboration 



2007 



Tyson] [2003) will introduce a huge leap in both sky coverage and 



in depth. These ambitious projects represent a culmination of the synoptic survey approach, and 
we will make a particular effort to examine their potential for supernova science. 



For our analysis, we will characterize each survey with four parameters 



1. the survey supernova depth, i.e., single exposure limiting magnitude mf^ for supernova de- 
tection when used in scan mode; this is set by collecting area (and monitoring time) 

2. the total survey scanning sky coverage, i.e., solid angle Ar^scan 

3. the scan revisit time ("cadence") Tvisit 

4. the total monitoring time Atobs? which (for a single cadence) is proportional to the total 
number Atobs /''"visit of visits 



For a fixed survey design and lifetime, these parameters are not independent, since exposure time 
comes at the expense of sky coverage and number of visits. 



There are numerous challenges and complexities in the process of e xtracting supernovae and 



their redshifts from surveys (and for sort i ng out their typ e s: see, e.g.. iDahlen fc Goobar 



Poznanski et al 



2007 



Tonrv et al.l (1200^ gives thorough discussion of these issue s wit h emphasis on Type la events; 



Kim &: Miguel 



2007 



Kunz et al. 



2006 



Blondin Sz Tonry 



2007 



Wan 



2002 



20071 ). 



Dahlen fcFranssonI (fl999l'). ICal-Yam fc Maozl (I2OO4I') and loda fc Totaiiil (|2005l V and the 



see also 



SDSS-II papers (jFrieman et al 



2008 



Sako et al 



20081 : Izheng et ahlbood l. 



Our simple survey parameterization cannot capture all of these subtleties, not do we intend it 
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to; rather, we wish our treatment to provide a rough illustration of the surveys' potential for core- 
collapse detection and science. Consequently, our parameter choices should be viewed as typical 
effective values, which may be different from (and weaker than) the raw survey specifications. 

For example, supernova identification and typing requires knowledge of the light curve. Thus, 
one c annot only ob serve the supernova at peak brightness, but also follow it after (and ideally be- 



fore). iTonrvl (|2003l ) recommends following the supernova for least 6m = l^^s below peak brightness; 
we will adopt this value as well. Thus, the effective supernova detection depth is mf^^^ = mmax — 6m, 
where mmax is the survey scan depth (i.e., depth for a single exposure). 

Note also that some upcoming surveys (such as DES) will only repeatedly scan a fraction of the 
sky which they map; but only the scanned regions will host the discovery of supernovae and other 
transients. Also, some surveys (e.g., Pan-STARRS and LSST) envision multiple periodicities and 
associated limiting magnitudes; for simplicity we will here chose conservative depths for the values 
given in Table [H to be consistent with the advertised scanning sky coverage. Thus one should bear 
in mind that in our analysis we have chosen the minimal parameterization one could use, which give 
only a simplified and idealized sketch of the real surveys. Given this, and the ongoing planning of 
future survey characteristics, our forecasts for the surveys' supernova results should be understood 
as indicative of the order of magnitude expected, but not as high-precision predictions. 



3. Core-Collapse Supernovae in a Cosmic Context 
3.1. The Cosmic Core-Collapse Supernova Rate: Expectations 

The total cosmic supernova rate (hereafter CSNR) 

^sn[^(W)] - J^'l (1) 

W •'com '^■'em 

is the number of events per comoving volume per unit time tem in the emission frame (i.e., cosmic 
time dilation effects in the observer's z = frame are not included). The total rate, and the various 
differential rates below, can of course be specialized to distinguish different groups of supernovae 
classified by intrinsic type and/or dependence on local environment. 

The present data on high-redshift core-collapse supernovae are too poor to construct an accu- 
rate CSNR. But th e CSNR is intimately related to cosmic star-formation rate = dM-i, / dVcomdt 



(jMadau et al.lll998l ). The connection is 



^SN = (2) 

where Xsn is the fraction, by mass, of stars which become supernovae, and (msN) is the average 
supernova progenitor mass (see Appendix [A]). A key point is that due to the short core-collapse 
progenitor lifetimes the two rates scale linearly, T^sN oc p^. The constant of proportionality depends 
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on the initial mass function (IMF). If the IMF changes with time (or environment) this comphcates 
the picture. In producing quantitative estimates we will follow most studies in assuming time- 
independent IMF. Thus the supernova/star- formation rate proportionality is a constant fixed for 
all time, namely TZsj^/p^ = 0.00915 (Appendix [A]) . 

Uncertainties in the cos mic rates for both s upernovae and star-formation remain considerable. 
As illustrated in detail by (jStrigari et al.ll2005l ). the cosmic star- formation rate is known to rise 
sharply towards redshift ^ ~ 1. In this low-to-moder a.te redshift regime, t h e shap e of the rate versus 
redshift is fairly well known, but as emphasized by iHopkins &: BeacomI (|2006l ) the normalization 
remains uncertain to within a factor ~ 2. At higher redshifts, the rate becomes even more uncertain, 
largely due to the paucity of data and also to uncertainties in our knowledge of the degree of dust 
obscuration. It is also worth noting that most studies to date directly or indirectly use massive 
stars as proxies for star formation. Consequently, the rate for cosmic massive star formation-and 
for cosmic supernovae-is less uncertain and IMF-dependent than the total rate. 

To illustrate the effects of these uncertainties on the synoptic survey supernova harvest, we have 
adopted two possible CSNR forms. These appear in Figuredl which shows the expected supernovae 
rate assuming a perfect environment (i.e. no dust extinction, etc). The s o lid cu rve in Figure H^a) 
is the C SNR derived from the cos mic star-formation rate of ICole et al.l ()200ll ) with parameters 
fitted by lHopkins &: BeacomI ([200a) (hereafter the "benchmark" CSNR). This rate sharply rises to 
a peak at z ~ 2.5, followed a strong but less rapid declines out to high redshift. To investigate 
the impact of the falloff from the pea k, we also show in the broken curve an alternate CSNR due 
to current observational data fitted by iBotticella et al.l (j2008l ) (hereafter the "alternative" CSNR). 
This rate also rise to redshift z ~ 0.5, though with a different slope; we somewhat arbitrarily set 
the alternative rate to a constant at z > 0.5 where the data are unclear; in any case we will find 
that few events from this high-redshift regime will be accessible to the all-sky surveys which are 
our focus. 

Synoptic surveys will measure several observables associated with cosmic supernovae: their 
numbers and location, and some portion of their light curves in different bands. Spectroscopic 
redshifts of host galaxies can also be determined (when visible; see §6.5p . Using the number counts 
and redshift indicators, one can deduce an observed core-collapse rate, per unit redshift and per 
unit time and solid angle. This observed rate distribution directly encodes the CSNR via 

(3) 



dNsN 



dNsN dtcm dVcc 



dQdtoiysdz dVcomdtcm dtohs dQdz I + z dz 

where Vcom is the comoving volume and rcomiz) is the comoving distance out to redshift z. The 
1 + z factor corrects for time dilation via dtobs = (1 + z)dti,m- 

Figure [IJb) shows the all-sky cumulative frequency of cosmic supernovae for an observer at 
z = 0, i.e., 



dNsN 

dt 



(< ^;)all-sky = 47r 



dNsN 

dildtohsdz 



dz' = 47r / 7^sN(z') 



drc 



1 + z' dz' 



dz' 



(4) 
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„ 0,0020 




redshift z 



Fig. 1. — (a) Top panel: Possible cosmic core-colla pse supern o va ra tes as a function of redshift. 

(I2OOII) cosrn ic star-formation rate; 



Cole et al 



The solid curve is the result calculated based on the 
the broken curve is based on current supernova data (jBotticella et al.l 120081 ): see Appendix [XI (b) 
Bottom panel: The idealized, all-sky cumulative rate of all supernovae observed over redshift to 
z, for an observer with no faintness limit and with no dust extinction anywhere along the line of 
sight. 
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These curves give the total rate of observed cosmic supernova explosions out to redshift z for 
an idealized observer monitoring the entire sky out to unlimited depth and without any dust 
obscuration anywhere along the line of sight. 

All of these idealizations will fail, some of them drastically, for real observational programs. 
Nevertheless, one cannot help but be tantalized by the enormous explosion frequencies indicated 
in Figure [Hb). With our benchmark CSNR, out to redshift z = 1, something like ~ 1 supernova 
explodes per second somewhere in the sky. Out to redshift z = 2, this rate increases to ~ 6 
events/sec. Clearly, even with a small detection efficiency, synoptic surveys are poised to discover 
core-collapse supernovae in numbers far exceeding all supernovae in recorded history to date. 

For numerical results in Figure [Hand throughout this paper, we adopt a flat cosmology with 
= 0.3 and Qa = 0.7. For the Hubble constant we adopt the value Hq = 71 kms~^ Mpc~^, i.e., 
h = 0.71 where Hq = 100 h kms~^ Mpc ~^. These values are consistent with recent determinations 



using WMAP and large-scale structure (jSpergel et al.l 120071 ) 



3.2. The Effect of Dust Obscuration 



The enormous inventory of cosmic supernovae is not, unfortunately, fully observable even 
for arbitrarily deep surveys. In a realistic environment there are several factors which will hide 
the supernovae from us; dust extinction is one of the most important, and probably the most 
uncertain. Core-collapse supernovae mostly explode within regions of vigorous star formation 
which are thus likely to be dusty environments. Consequently, we expect that some core-collapse 
supernovae will be obscured to the point where they are not detected in synoptic surveys. The 
fraction of supernovae lost to dust obscuration, and particularly the possible redshift dependence 
of this extinction, represents a crucial systematic error which must be addressed before one can use 
survey data to infer information about supernova populations and their cosmic rates. 

For the purposes of our present estimates of survey supernova yields, we follow the approach of 
Mannucci et al.l (120071 ). These authors characterize losses due to dust extinction and/or reddening 
in the host galaxies via a fraction adust (-2) of undetected events at each redshift. This fraction could 
in principle differ for the various core-collapse types; for the present treatment we will assume it is 
the same for all such events. As core-collapse statistics become available from surveys, this issue 
can and should be revisited; more on this below, and in ^ The resulting fraction of detected 
supernovae is thus the complement /dust (-2) = 1 — adust(^)) which measures the reduced supernova 
detection efficiency in the presence of dust. Expressed as an effective extinction A for the supernova 
population at z, we have Aes{z) = —1.086 In /dust- 



Mannucci et al.l (120071 ) estimate the fraction of missing supernovae by comparing the observed 



detections in the optical with those in radio and near-IR. They conclude that dust evolution is very 
strong; this becomes a dominant limit a .tion t o the discovery of core-collapse events at high redshift. 
In the local universe, iMannucci et al.l (|2007l ) find that the vast majority of the events occurring in 
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massive starbursts (luminous infrared galaxies) are missed. Because these galaxies harbor only a 
small fraction of the local supernova population, the overall optically missing fraction at z = is 
estimated to be rather modest: Odust = 5 — 10% . If, however, high-redshift star formation occurs 
in starb urst environrnents (i.e., luminous and ultra-luminous galaxies , which are highly exti ncted; 
see e.g. 



Choi et al 



Small et al.1 Il997l : iHuehes et al.1 Il998l : IPerez-Gonzalez et al.1 l2003l : iLe Floc'h et al.1 12005 : 



20061 ) then the fraction of missing events rises sharply with redshift. Multiwavelength 



observations of light from pr e-supernova massive stars a lso supports the idea of increasing dust 
obscuration at high redshift. lAdelberger Steidell (j2000l ) find that ultraviolet light from massive 
stars in z ~ 3 galaxies is mostly reprocessed by dust into thermal submillimeter emission, so that 
the observable galaxy luminosities have Lsub~mm/-^uv ~ 1 — 100. 



Mannucci et al.l (|2007l ) estimate the portion of supernovae which will be "catastrophic losses" 



to severe extinction, and propose that the missing fraction can be described by a linear relation 
adust(-2) = 0.05 + 0.28z for the core-collapse supernovae for redshift z<2. Thus, the fraction of the 
supernovae which remain optically d etect able is fdnst(z) = 1 — adust(-z) = 0.95 — O.282; for z < 2. 
At higher redshift. IChen et al.l (120071 ) and lGnedin et al.l (j2008l ) argue that /dust is small; they find 
limits consistent with /dust = 0.02 for these redshifts. 

We will smoothly match these two estimates, and adopt a fraction of the supernovae which 
can be detected after dust extinction of 



/dust(^;) 



0.95 - 0.28z , z < 3.3 
0.02 , z > 3.3 



(5) 



For these values of /dust) the effective extinction varies from ^eff = 0.056™**^ at z = to A^s = 
^_25™ag 2 > 3.3. In practice, we will find that cosmological dimming of supernovae beyond z ^ 1 
is itself so large that surveys up to and including LSST will see relatively few events, so that the 
details of the adopted dust model in this regime will not affect our conclusions. 



The strong redshift evolution of dust obscuration in the empirical iMannucci et al.l (120071 ) model 
deserves comment. Prom a physical point of view, the rise in dust losses a(z) towards high redshift 
implies that at earlier times, the birth environments of supernovae are significantly more enshrouded 
than those now. This interesting result itself deserves a deeper elucidation, one which will likely be 
easier to formulate and test in the presence of survey supernova data. Prom more practical point 
of view, our adoption of a model wherein dust losses grow rapidly with z should yield conservative 
(or at least not optimistic) predictions for the supernova harvest at large redshifts. That is, if it 
turns out that host galaxy effects do not change rapidly with cosmic time so that the efficiency of 
supernova detection remains close to the high local value, then our rate predictions at z ~ 1 would 
be boosted by a factor of ~ 1.5. 



Note also that /dust as IMannucci et al.l (|2007l ) and we have defined it characterizes the observ- 
able portion of the ensemble of supernovae at a particular redshift. Implicitly, individual supernovae 
are treated as either detectable or not, i.e., dust effects are considered negligible or total; our cal- 
culation treats total, catastrophic losses of supernovae using this /dust formalism. We separately 
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include the effect of partial extinction due to dust, where the apparent magnitude of supernova 
is reduced but still visible, as discussed in the next section. Of course in reality, all supernovae 
will experience some level of extinction in their host galaxies, with the distribution of host-galaxy 
extinctions changing with redshift. A more detailed study of dust effects on supernovae (and uses 
of supernovae to quantify and calibrate these effects) would be of interest for further investigation; 
see discussion in ^ 



3.3. Supernova Observability at Cosmic Distances 

3.3.1. The Supernova Luminosity Function 

Locally observations of core-collapse supernovae reveal diverse light curves, with a wide range 
in peak luminosity, and very different evolution after maximum brightness. The vast majority 
of supernovae discovered by synoptic surveys will lie at cosmological distances, and thus will be 
detectable mostly near their maximum luminosity. Thus we will focus on the observed distribution 
of peak brightness, and the timescales on which supernovae sustain it. 

The distribution of peak absolute magnitude Mpeak is given by the supernova luminosity func- 
tion (/)snif,x = 0snif,x(-^peak; z) which may have a redshift dependence; we choose a normalization 
such that at any z, J i;^snif,x(-^i dM = 1, with this, we may write the cosmic supernova rate per 
absolute peak magnitude as 

— ^f^^^ = 7^SN(^) </'snlf,x(Mpeak) (6) 

WVcom "■(-em "-iKipeak 

where here and throughout the possible redshift dependence of the luminosity function is under- 
stood. 



Richardson et al.l (|2002l ) find the best-fit formulae for the supernova peak luminosity functions 
in i5-band for different types of supernovae based on their tabulation of 279 supernovae of all 
types, for which absolute magnitudes were available at peak brigh tness. Of these, there w ere 168 



events of all core-collapse types: II-P,L,n and I-ab. For each type, iRichardson et al.l (|2002l ) fit the 
observed S-band absolute magnitude distributions with gaussian profiles, in some cases including 
two profiles where the data suggested "bright" and "dim" subclasses. Their results provide the 
basis for the luminosity functions used in this paper. 

Note that we use the observed distributions rather than intrinsic, dust-corrected versions. Thus 
we automatically include the mean extinctions (ranging from A ~ 0.1™'^^— O.S™'^^ for different types) 
found at low redshift. The redding effect due to dust (ranging from E ~ 0.02 ™'^^ — 0.26™'^^ acros s 



different bands and redshifts) is also added, based on the information given bv lKim Sz Led (|2007l ). 
As noted in the previous section, catastrophic losses of supernovae due to large extinction and its 
possible evolution at high redshift is treated separately via our /dust parameter. 
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We adjust the iRichardson et alj (|2002l ) distributions in two ways. First, we converted from 
their Hubble constant of h = 0.6 ot our adopted value h = 0.71. More importantly, we assume 
that each gaussian is a good representation of the data around the peak, but we do not allow 
the wings to extend arbitr arily far. Instead , we cu t off the distributions at \M — MmeanI > 2.5cr, 
where no data exist in the IRichardson et al.l (|2002l ) sample. We introduce these cutoffs in order to 
avoid extrapolating to very rare, bright events which in a large survey could extend the redshift 
reach considerably. Below ( §4.2p we discuss the effect of this cutoff and its effect on the predicted 
supernova redshift range. 



3.4. Supernova Discovery in Magnitude-Limited Surveys 



Surveys will discover supernovae monitor lightcurves in one or mo re passbands Here we will 
adopt the SDSS ugriz photometric system, which uses AB magnitudes (jFukugita et al.lll996l ). 



The light curve of any supernova will suffer redshifting and time dilation effects. For passband 
X we put 



= Slog 



V do 



+ K^{z) + A^iz) = ii{z) + K^{z) - 1.086 In /< 



dust 



(7) 



with di the luminosity distance and ^{z) is the usual distance modulus with do = 10 pc. The dust 
extinction A is included via the factor /dust (eq. E]). The X-correction accounts for redshifting of 
the supernova spectrum, and is discussed in Appendix ICl 

As noted above, at each redshift the effect of dust will be to obscure some fraction of supernovae. 
The remaining unobscured events will have apparent x-filter magnitudes oirux = Mx+ii{z)+Kx{z). 
The expe cted m.-^ distributi o n thu s reflects the underlying distribution of absolute magnitudes Mx ■ 
Since the IRichardson et al.l (|2002l ) supernova luminosity function we use is in the -B-band, we need 
to find the corresponding S-band magnitude in order to find the right corresponding number of 
supernovae; this transformation to is straightforward and is given by rux = ms + tJxB, where 



J^f F{\)Sx{X)dX 

rjxB = —2.5 log — ^ h zeropomt correction 

/|/F(A)5B(A)dA 



(8) 



is a color index which translates between the x and B magnitudes in the rest frame, and zeropoint 
correction is the correction for different zeropoint of the SDSS magnitude syst em and the Johnson 
magn itude system. For the spectral shapes F(X) we use the prescriptions of iDahlen Fransson 
(|l999l ) as described in Appendix [Q 



The absolute x-band magnitude distribution of unobscured supernova e at redshift z is (t><in}f,^\M x- 
Vxb], where i;^snif,x is the luminosity function in i?-band as tabulated by IRichardson et al.l (j2002l ). 
Therefore the distribution of a certain type of supernova apparent magnitudes rrix in x-filter is 
0snif ,x ["T-x — fJ-iz) — Kx{z) — rjxB], Thus the fraction of all (unobscured) supernovae at z which 
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fall within the survey x-h&nd magnitude limit is a sum over the luminosity functions for all 

core-collapse types: 

, ( \ _ ST- /""""^ (psnKA'^x - Kz) - Kx{z) - r]xB\ dm 

/maglimi^j - 2^ ( dm. 



types 



E 



/ 4'snU,x{rn) dm 



types / <^snlf,x("i) dm 

= /snif[< Mu,n{z,ml^^)] 
which is the cumulative fraction of supernovae whose absolute magnitude is brighter than 

Mn^{z, mf^J = mfi^ - /i(z) - K^{z) - (9) 

To develop some intuition, suppose the supernova peak brightnesses lie in a range Mpeak £ 
(Mbright) -^dim); and ignore for now the effects of dust. Then for low redshifts such that the absolute 
magnitude limit Miim from eq. Q is fainter than Mdimi we can expect to see all supernovae, and 
/magiim = 1- For these redshifts, we can study the entire supernova luminosity function and test 
whether it varies with redshift. On the opposite extreme, for high redshifts such that is dimmer 
than Mbright we can see no supernovae, so /magUm = 0; this then defines the survey redshift cutoff (for 
fixed mf^). Finally, for intermediate z such that mjj^— Mjim < fjL{z)+Kx{z)+r]xB < Mbright! 
we have < /magUm < 1 ; over these redshifts the survey samples the bright end of the supernova 
luminosity function. 

Both magnitude limit and dust extinction reduce the expected supernova detection, and do so 
independently of each other. Consequently, we can find the net supernova detection probability by 
simply taking the product of the individual factors: 

/detect (^j '^li^n) = /magUm(^j ^^fi^n) /dust(^) (10) 

Figure [2] shows the resulting detectable fraction of supernovae. The left panel shows the shape 
of /magUm for the g and r bands. At redshifts close to zero, /magUm ~ 1 which means that almost 
all supernovae are detected in the local universe. And it approaching to zero at high redshift, 
which reflects the fact that no supernovae can be detected at high redshift because of the survey 
deepness. Note that g and r bands are competitive for mf^ < 24, but for higher mf^, 
g-h&nd drops a lot faster than those in r-band especially around z ~ 0.4, which is cause by the 
effect of K-correction. The figure also shows that for higher mf^, /magUm decays less rapidly. The 
right panel shows /detect for different mf^^, using our adopted dust model (eq. [5]). We see /detect 
shows the same trend as /magUm except the detectable fraction is reduced due to dust and we can no 
longer observe all supernovae even in the local universe. It is also clear to see that going to fainter 
mjj^ significantly boosts the detectable fraction at high redshift. For mjj^^ = 23™*^^, /detectable is 
almost zero at redshift z ~ 1 for both g and r bands. But going to m^^ = 26™'^^, ~ 55% of the 
supernovae at redshift z ~ 1 remain visible both the g and r bands. 
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0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 

redshlft z redshlft z 



Fig. 2. — (a) The fraction of supernovae detected based on different survey deepness in g and r 
bands, with ranging from 21™^^ down to 27™^^; effects of dust obscuration are not included, 
(b) As in (a), but including the effects of dust obscuration strongly evolving with redshlft as modeled 
by eq. ([5]). 

This means that deeper surveys (and/or scanning modes in which smaller areas are scanned 
more deeply) will probe supernovae out to much higher redshifts. Deep survey modes will also 
probe a much wider regime of the supernova luminosity function and light curves over a broad 
range of cosmic epochs, thus testing for redshift evolution in supernova properties. The clear lesson 
is that the scan mf^ is critical in determining the quality and reach of the supernova science. In 
particular, we urge that scans strategies include modes which push > l^'^s deeper than the all-sky 
depth. 



3.4-1- Supernova Light Curves 



The observed population of core-collapse superno vae shows a broad range of timescales and 
time histories in their decline from peak brightness (e.g., Doggett &: Branchlll985 : Leibundgut Sz Suntzefl 
20031 ). The amplitude and time behavior of these curves encodes a wealth of information about 
the underlying physics of the supernovae as well as their interaction with the circumstellar and 
interstellar medium. 

Empirically, light curves broadly fall into phenomenological categories, those whose magnitudes 
decline in a relatively steep, linear way (Type II-L ) and those which linger near peak brightness 
with a relative plateau in magnitu d e (Ty pe II-P). iPatat et al.l (|l993l ) compiled 51 Type II light 
curves, and analysis in iPatat et al.l (| 19941 ) showed that plateau-type supernovae typically decline 
from peak brightness at rates which vary the range (O.?™*^^ — S.l^'^sjy'ioo days, while linear-type 
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events typically have {3.9^^^ — 5.7™'^^)/100 days. Unfortunately, the lightcurves available at the 
time of these studies were poorly sampled near the peak itself, where the behavior is most critical 
for our purposes. 

Fortunately, subsequent data, particularly using Swift, give s a clearer picture of th e early light 
curves for a few events. For plateau event SN 2005cs, data in iPastorello et al.l ([200a) show that 
~ 15 days after peak brightness, the supernova dimming was strongly depending on passband: 
AMi5(C7) ~ l.S^'^s, AMi5(B) ~ 0.7"^^s, AMi5(y) ~ O.lS^'^s, and AMi^iR) ~ O.l'^^s. Another 
Type II-P event, SN 2006bp, a fter ^ 13 days dec lined by ~ in U, but within errors was 

essentially constant in B a nd V (iDessart et al.ll2008l ). For Type lb, the recent event SN 2008D was 
seen from shock breakout (|Modiaz et al.ll2008l ): after dropping from this brief initial outburst, the 
flux increased for ~ 15 days to a maximum. Afterwards, the brightness decline rates lengthen with 
wavelength, with a drop of AM ~ after ~ 10 days in [/-band, but after about 15 and 20 days 
in B and V respectively. 

These multicolor data show that brightness decline in V and longer passbands comparable t o 
if not slower than the typical range of AM15 ~ l"^^s _ 2™^s found in Type la events (jPhillipslll993l ). 
This implies that surveys timed for Type la discovery will automatically be well-suited and possibly 
even better-sampled for core-collapse events. In particular, we will find below that the r and also 
g passbands are the most promising for survey supernova detections. Thus, if cosmic supernovae 
follow the behavior of these local events, we expect that the light curves will remain within, e.g., 
Am ~ 0.5^*^^ of peak brightness (a factor 1.5 in flux) for a timescale of at least a week. In some 
cases this timescale will be longer, and possibly also with detections in the rising phase. 

For synoptic surveys to detect core-collapse supernovae near their peak brightness, the cadence 
needs to be shorter than the (observer-frame) brightness decline time. Thus weekly revisits are 
sufficient for marginal detections, and cadences of ~ 3 — 4 days will often see the event three or 
more times. In the cases of plateau events, the supernova should remain near peak brightness 
for many such revisit times. Furthermore, due to cosmological time dilation effects, the observed 
brightness decline timescale Tobs = (1 + •z)''"rest is increased by a factor of 1 + z, which extends the 
detection window and offers a greater opportunity to recover a well-sampled lightcurve. Also, we 
see that color evolution is not strong in V and R bands. The UV and blue do fade more rapidly, 
and the supernova reddening depends on the type. For events where bluer rest-frame colors are 
available, this might be a useful means of photometrically determining supernova type. 



4. The Cosmic Core-Collapse Supernova Rate: Forecasts for Synoptic Surveys 



In this section we will work out general formalism for supernova observations by synoptic 
surveys. We then apply this formalism to specific current and proposed surveys 
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4.1. Connecting Cosmic Supernovae and Survey Observables 

4-1.1. General Formalism 

It is useful to define a differential supernova detection rate per unit redshift, solid angle, and 
apparent magnitude in x-band: 

dA^SN,obs,x „ , . r{zf dr l \ t \ l 

- TZsN{zj — -T- Jdust{z) (psnlf,xl?7lx - - Kx{z) - Tj^B] (H) 



dtobs dz dQ dm 1 + z dz 

This expression adds the effects of supernova luminosity (cf eq. [6]) and of dust obscuration (eq. [5]) 
to the ideal rate of eq. Throughout, we will for simplicity refer to the entire core-collapse 
supernova population, but the formalism could equally well distinguish the va rious core-collapse 



types, and compute the rates of each. An example of such a treatment is the IScannapieco et al 



(|2005l ) study of the rate and detectability of pair-instability supernovae. 



The differential rate in eq. ()lip relates the observables in a synoptic survey to underlying 
properties of cosmic supernovae. As such, a wealth of information can be recovered by a good 
statistical sample of supernovae over a redshift range: one probe different terms and their underlying 
physics. For example, at fixed z, the range of observed supernova magnitudes in x-band probes 
the supernova peak luminosity function (/)snif,x(-^; z) at magnitudes Mx = rux — fi{z) — Kx{z) — rjxB- 
Comparing these results at different redshifts with local determinations can reveal any redshift- 
and/or environment-dependence in the core-collapse supernova luminosity function. 

Another aspect of cosmic supernovae probed by synoptic surveys, and central focus of this 
paper, is the cosmic supernova rate. Whereas the supernova luminosity function can be determined 
from the distribution of supernova magnitudes at the same redshift, the cosmic supernova rate 
comes from the distribution of supernova counts across different redshifts. The observed differential 
rate for supernovae of all magnitudes in the x-hand is 

rsN,obs,xU; = -J. = / dm — = 7^sN(^) /detect ,x ; — — — (12) 

otobs dzdil J dtohs dz as I dm 1 + z dz 



Note that this is the idealized rate of eq. ([3]) reduced by the detection in x-band /< 



detect ,x- 



One can get a sense of the orders of magnitude in play via the definition of a dimensionful 
scale factor 

TsN.o = ^sn(O) 4 = 7.2 X 10^ events yr'^ sr^^ -s i (13) 



10-4 yr-1 Mpc- 

0.22 events sec- (14) 
2.2 X 10= events y.-'deg-(^^,,^pf(|-,) (15) 



We may then define a dimensionless distance u{z) = r{z)/dH, with dn = c/Hq the Hubble length. 
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Fig. 3. — The cosmic supernova detection rate in r-band, expressed in number of events per solid 
angle per time, shown a function of redshift. The curve labeled "unobscured" ignores both effects 
of dust extinction or the flux limit of the survey (i.e., /detect = !)• The curve labeled "with dust" 
includes dust extinction only, but with rrv^ = 00. The remaining curves are for surveys with rrV]^ 
as labeled, and include dust extinction. Note that the vertical axis is shown both in units of events 
per second per steradian (left scale) and events per year per square degrees (right scale). 



and write 



/ ^ _ -r ^sn(^;) u{zf du 

S'N,obs\Z) — i SN,0 7^ Z — ; /detect, 

7v.sn(0) 1 + z dz 



sn \ 
lim/ 



(16) 



Figure [3] plots the observed supernova rate FsN.obs per solid angle in r-band. For comparison, 
we show the idealized cases of mf^^ = 00 and /dust = 0, as well as realistic cases in the presence of 
dust and with different mf^. The amplitudes of the curves in Figure [3] confirm the large numbers 



of events expected from eq. (fT3|) . 

The shapes of the curves can also be readily understood. At low redshifts, the surveys see most 
of the supernovae that occur-i.e., the entire luminosity function is sampled; cf Figured! Hence at 
small z, the supernova sample is simply limited by the cosmic volume within z: F oc dVcom/dz ~ 
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ftoradfcoia/ dz ~ Thus the detection rate initially rises quadratically with z; this volume effect 
is essentially independent of survey magnitude limit, as we see by the overlap of the curves in this 
regime. 

In the high redshift limit, several effects act to suppress supernova detectability. At z > 1, Tcom 
rapidly saturates at the comoving horizon scale, and nearly all observable cosmic volume is sampled; 
in this regime, the volume factor decreases as dVcom/dz ~ drcom/dz ~ 1/H{z) ~ (1 + z)"^/^. In 
addition, time dilation effects become large and add another factor of (1 + 2;)^^. For these reasons, 
even the idealize d (unobscured, Ta\f^ = 00) rate drops. Moreover, in some models (such as that of 



Cole et al.ll200ll ). the CSNR itself is intrinsically expected to drop after a peak, perhaps somewhere 
in the range z ~ 1 — 3. On top of this, the effects of dust obscuration become large at z > 1 and 
removes further supernovae in this range. Finally, a finite survey magnitude limit truncates still 
more events at high z. 

The combination of the low-redshift rise and high-redshift drop acts to create a peak in su- 
pernova detectability. The position of the peak is sensitive to the CSNR itself, and the details of 
dust obscuration. But the peak position and amplitude are also both very sensitive to the survey 
magnitude limit; both rise sharply as survey depth rrV]^ increases. This illustrates a key conclusion 
which will be manifest in several other ways below: for discovery of core-collapse supernovae at 
high redshifts, the most important aspect of a synoptic survey is its limiting magnitude; investment 
in deep scan modes (m^^ > 24 mag) will reap substantial rewards. 

Figure Hlshows the same supernova rate redshift distribution as in Fig.[3l but for the five ugriz 
passbands with SDSS filters and efficiencies. For each band we fix m^?^ = 24™'^^, gge that the 
discovery rate is the highest in r for essentially all redshifts, with 5-band counts very nearly the 
same except around the peak at 0.2 < z < 0.6. The relative smallness of the counts in other bands 
traces back predominantly to low detector efficiency in i and z, and redshifting effects for u. The 
upshot is that for synoptic surveys, r and g bands are (in that order) clearly the most promising 
for supernova search. 

We have thus far shown the total supernova rate redshift distribution, summed over all core- 
collapse subtypes. Figure [5] illustrates how the different subtypes contribute to the aggregate. Here 
we fix Tn^f^ = 24^*^^ and show results for t he r and g bands. It is worth recalling that we have 



assumed the low-redshift iRichardson et al.l (|2002l ) determination of luminosity functions and type 



distributions holds for all redshifts. In this scenario, we see that in both bands. Type Iln events 
give the largest contribution to the signal at z > 0.3, and totally dominate the counts at z > 0.6. 
This is expected, since it is the intrinsically brightest core-collapse subtype. Thus the redshift reach 
of supernova discovery (and associated results such as the CSNR) in synoptic surveys will depend 
sensitively on nature Type Iln events at z > 0.6. It will thus be crucial to determine whether 
these events show evolution in their luminosity function and/or relative fractio n of core-collapse 



events with redshift (e.g., via metallicity effects). Also, it is worth noting that the lRichardsoii et al 



(|2002l ) luminosity function we have used is relatively narrow. As noted recently bv ICookd (|2008l ) 
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Fig. 4. — Number of supernovae per year per solid angle per redshift with = 24™*^^ in different 
bands. 



some Type Iln events have now been observed with luminosities far above the range of values we 
consider. If so, then the redshift range of synoptic surveys could thus extend significantly further 
than in our estimates. 

Figure [5] further predicts that the other core-collapse types should have observably distinct 
redshift ranges in different bands, again assuming no evolution in luminosity function or type 
distribution. The upper panels of Fig. [5] show the individual subtype detection rates, as well as 
their sum. Type II-L events have simlilar behavior in both r and g bands, peaking at z ~ 0.45 
then rapidly dropping o ff. Although Type II-P events are the largest core-collapse subtype in the 
Richardson et all ( 2002 ) sample, they are also by far the intrinsically dimmest, ~ i^^g _ gmag 
fainter than the other types. We thus find that Type II-P have a smaller redshift range than Type 
II-L and Iln events. The counts and redshift range of Type lb and Ic events are notably different 
in the two passbands. This traces to the effects of UV lineblanketing which removes blue flux; thus 
at high redshift the i^-correction first shifts photons out of the g band, with the r-band signal 
surviving until higher redshift. Note also that the "bright" and "normal" Type lb and Ic events 
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Fig. 5. — Supernova rate redshift distribution, as in Fig. [H broken down by core-collapse type. 
Results shown for (a) r band, and (b) g band; both have mf^^ = 24^*^^. Top panels: detection rate 
distribution per subtype; bottom panels: fraction of each subtype rate relative to the total. We see 
that intrinsically bright Type Iln events dominate the counts at high redshift (z > 0.5) and thus 
determine the redshift reach for core-collapse discovery. 



lead to the double-peaked structure in their redshift distribution. 



The lower panels of Fig. [5] shows our forecast for subtype fraction detected as a function of 
redshift, i.e., the ratio of each subtype rate to the total. At z = 0, t he subtype fractions go to 
the observed local values we have adopted from iRichardson et al.l (|2002l ). as required by our model 
design. For z ~ — 0.2, we see that all subtypes make significant contributions tot the total, 
and thus for this redshift range, the sharp rise in the total detection rate (top panel) is due to 
contributions from all subtypes. The features around the maximum in the total rate {z ~ 0.2 — 0.5) 
are due to the interplay between the rise of the Type Iln events and the successive dropout of the 
other types. Finally, we see that for z ~ 0.5, Type Iln events essentially completely set the total 
rate. 

Because the highest-redshift detections will be dominated by Type Iln events, the nature of 
and evolution of this subtype will play a cruci al role i n sett ing the high-redshift impact of surveys 
for core-collapse events, as also pointed out bv ICookd (j2008l ). As we have noted, intrinsic evolution 
of the Type Iln fraction of core-collapse events would directly change-and be written into-the high- 
redshift signal. But at present, the uncertainties are very large even when evolution issues are set 
aside. Namely, published data are as yet very uncertain concernin g the local, z ^ fraction of core- 



collapse events which explode as Type Iln. Our forecasts use the IRichardson et al.l (|2002l ) sample 
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which finds 9 Type Iln events out of 72 core-collapse events, for a fracti on of 12.5%. Howev e r, thi s 



discovery fraction are very uncertain. For example, the prior work of iDahlen &: FranssonI (|l999l ) 
compiled their ow n core-collapse discovery statistics, and adopted a Type Tin event fraction of 2%, 
while noting that ICappellaro et al.l (|l997l ) recommend a Type Iln fraction of ~ 2 — 5%. Because 
the high-redshift core-collapse detections will be dominated by Type Iln events, if these values 
better reflect the intrinsic fraction, this would dramatically reduce our predicted detection rates for 
z > 0.5 by factors of ~ 2 — 6, and thus also reduce the maximum redshift at which core-collapse 
events can be seen in surveys. Clearly, the small numbers available when all of these compilations 
were made render the Type Iln fraction estimates uncertain; indeed, to a lesser extent the estimates 
for the more common core-collapse types suffer similar problems. 



In light of the uncertainties in the iRichardson et al.l (j2002l ) and prior compilations, it is worth 
noting that considerably more supernova data already exists. A detailed, systematic study of the 
luminosity function and intrinsic subtype fractions of local events would be of the utmost value for 
forecasts of the sort we have presented. Moreover, precise and accurate local measurements will 
play an essential role as a basis of comparison for the future medium- to high-redshift data, in order 
to empirically probe for evolution within and among the core-collapse subtypes. 



4-1-2. Unveiling the Cosmic Core-Collapse Supernova Rates 

As noted above, synoptic surveys will revolutionize our understanding of the CSNR because 
they will directly determine the rate through counting- We now are in a position to determine 
the supernova counts for realistic (magnitude-limited, dust-obscured) surveys. Using these, we can 
demonstrate how the CSNR can be extracted. We can further determine its statistical uncertainty 
and the impact of survey depth and sky coverage. 

Consider a survey with scan area AOscan and limiting magnitude ^ij^, the total number of 
supernovae seen in x-band in time Atobsi in a small redshift bins of width Az = zj — Zj <C 1 centered 
around z = {zj + Zi)/2 is 

AA''sN,obs,x = A$7scanAtobs rsN,obs,x(2:) (17) 

A A '^Sn(-2) u{zY du t sn \ 

= AiZscanAtobs Az TsN.Cx ^^^^^-^ I + z 'dz J'^^^^^^^^^^'''^^^^' ^ ' 

Thus we see that the cosmic supernova rate is directly encoded in our binned data. This means we 
can use the binned data to extract the supernova rate: 

,obs,x /irv\ 

''^^^^) = Aa;:: Ai^d^^---^^^"^-) 4A^ 

this result is a major goal of this paper. Physically, we see that as we accumulate supernovae, i.e., 
as AAsN,obs fills out the redshift range accessible to the survey, we obtain an ever better measure 
of the SN rate. 
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We can also compute the statistical uncertainty in the CSNR derived from counts in surveys. 
The statistical error arises from the counting statistics in the supernova number. Expressing this 
as a fractional error, we have 

CtCR-sn) _ cr(AiVsN,obs,x) _ 1 _ ^20) 



obs,x 



But from eq. (|18p . we see that AA'^gN obs,x scales linearly with the product of detected fraction and 
survey sky coverage, as well as monitoring time and redshift bin width. Thus we find the CSNR 
statistical error should scale as 

^(^sn) _ 1 ^21) 



scan Atobs Az 



V /dctcct.xV^^j '^liJn) Atobs AJlgcan 

Consequently, for a fixed redshift bin size Az, the CSNR accuracy grows with the product Aiobs Af^scE 
and implicitly with via the detection fraction. Thus survey sky coverage and magnitude limit 
(i.e., collecting area) enter together, and we see the payoff of a large survey etendue. 

Thus, we can find the survey properties, needed to achieve any desired precision in the CSNR 
at some redshift z. For a fixed nif^ and thus /detect; monitoring time and sky coverage enter 
together as the product Atobs AQgcan- Figure [6] shows the needed monitoring time Atobs Afigcan to 
measure the CSNR to a statistical precisions of (Tstat('^SN)/^SN = 10%, and with different survey 
mf^ in r-band. In both panels we choose Az = 0.1 for the redshift bin size. The two panels show 
our baseline and alternative CSFR. From these figures we can see that these two different adopted 
CSNR behaviors both yield very similar results for the survey CSNR detectability. 



Again the shapes of the curves can be understood. As shown in eq. (j2ip . that the precision at 
each bin scales inversely with the supernova differential redshift distribution as robs(-2)~^''^- Not 
surprisingly therefore, the least monitoring is needed to measure the CSNR for z near the peak in 
the redshift distribution On the other hand, redshifts in the high- and low-redshift tails of rgN,obs 
require increasing monitoring, eventually to the point of unfeasibility. 

Figure [6] makes clear that increasing mf^ brings a huge payoff reducing the needed monitoring 
At Ar^scan- To achieve a (t(7?.sn)/'7^sn < 10% precision at redshift z = 1, the monitoring becomes 
about 1000 times smaller in r-band if we increase from 23™'^^ to 26™^^. Clearly, for any survey, 
increasing will drastically shorten the observing time needed for the high redshift supernovae. 
In practice, given fixed survey lifetimes, this means that sets the maximum redshift reach over 
which the survey may determine the CSNR (via eq. [9j). 
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4.2. Forecasts for Synoptic Surveys 



For a given survey witli a fixed scanning sky coverage AOscan; we can determine the total 
number of supernovae expected in each redshift bin. We can also forecast the accuracy of the 
resulting survey determination of the CSNR. Namely, we can turn our sky coverage-monitoring 
time result (Figure [6]) into a specific prediction for the needed time to determine the CSNR to 
a given precision. In practice, this amounts to a determination of the redshift range over which 
different surveys can measure the CSNR. Our detailed predictions appear in Appendix (Bj here we 
summarize the results. 
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Fig. 6. — Survey CSNR discovery parameter, i.e., the product of survey monitoring time and sky 
coverage At Ai} needed to measure the CSNR to a specified precision. Data are binned in redshift 
units of Az = 0.1 vs. redshift. Top panel: discovery parameter needed to reach 10 % precision with 
our benchmark CSNR. Bottom panel: discovery parameter needed to reach 10 % precision with the 
alternative CSNR seen in Fig. [TJ 
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Several main lessons emerge from considerations of specific surveys. When Pan-STARRS and 
LSST are online, these surveys will collect a core-collapse supernova harvest far larger than the 
current set of events ever reported. This alone will make synoptic surveys a transformational point 
in the study of supernovae. 

Moreover, synoptic surveys will detect core-collapse events over a wide redshift ranges. Table 
[2] summarizes the supernova redshift ranges correspond to the most likely mf^ of the surveys. The 
total supernova harvest depends sensitively on the survey depth, and in Appendix[B]the sensitivity 
to is shown. To determine the redshift ranges shown in Table [2l we set an (arbitrary) lower 
limit on the number of total supernova counts at A^min — 10* ^^e choose a lower redshift limit ^min 
such the cumulative survey supernova count in one year is A'survey(< -^min) = A^min- Similarly, the 
upper limit ^max is set by A^(> Zmax) = A^min is the number of supernovae detected within redshift 
z = Zmin within a year. 

As seen in Tables [2] and [3l the future surveys will find abundant supernovae over a wide 
redshift range. At low redshifts, the surveys will detect nearly all of the supernovae within the 
nearby cosmic volume accessible in their sky coverage. So surveys with a large AOscan) such as 
LSST and Pan-STARRS, have 

Zmin which docs not depends on itI]^- For DES, Afiscan is not 
as large, so that the number of supernovae brighter than rrV]^ = 21™'^^ does not accumulate to 
A^(< •2min)=10 Until Zmin=0.081. But for depths fainter than m\f^ = 24^**^, the survey does become 
volume-limited and the supernova counts accumulate to 10 at the same redshift. The upper limit 
of the redshift Zmax depends not only on sky coverage but also survey depth. For planned survey 
depths, DES will gather core-collapse supernovae to about z ~ 1.20; LSST will extend to z ~ 0.89, 
and could go further in modes with smaller sky coverage but deeper exposure. 

The large supernova counts and wide redshift ranges together mean that surveys will, by 
direct counting, map out the CSNR to high precision out to high redshifts. Future surveys should 
easily achieve 10% statistical precision for the CSNR for redshifts around which the survey's counts 
peak. We see that, with ra\^^ = 23™**^, LSST will reach out to z ~ 0.89, presuming that the 
relative fraction of the brightest, farthest-reaching events (of Type Iln and Ic) do not evolve with 
redshift. If so, then by direct counting future survey should witness the sharp CSNR rise. With 
deeper exposures and corresponding increases in redshift reach, surveys could begin to test for the 

Table 2: Survey Discovery Potential for Core-Collapse Supernovae in r-band 



Survey 


Expected Total 1-year 


SNII Redshift 


Name 


SNII Detections 


Range 


SDSSTI* 


1.70 X 10^ 


0.03 <z< 0.37 


DES 


2.74 X 10^ 


0.06 <z< 1.20 


Pan-STARRS 


5.14 X 10^ 


0.01 <z< 0.89 


LSST 


3.43 X 10^ 


0.01 <z< 0.89 



Note: *Reflects SDSS-II supernova scan season of 3 months per calendar year. 
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behavior of the CSNR above z = 1, a regime that is currently poorly understood. 

Both the survey yields of supernova discoveries, as well as their redshift ranges, are strong 
functions of survey depth. As shown in Table [3] of Appendix[Bl each magnitude increase in survey 
depth yields a large enhancement (a factor ~ 3) in total supernova counts. This in turn leads to 
large enhancements in redshift range, and thus in the range over which the CSNR is measured. As 
shown in Appendix [Bl increased monitoring time needed to achieve higher mf^ will come at some 
cost, though this will be partially offset by the higher supernova yield in a deeper exposure. Finally, 
for the large population of low-redshift supernovae, deeper surveys will lead to better lightcurve 
determination, allow for a more accurate photometry over a larger brightness range and thus longer 
timescales. 

As noted in §3.3.11 our fiducial results are for supernova peak magnitudes whose luminosity 
functions (each of which is one or two gaussians for each core-collapse type) are nonzero only 
within \M — Mmeanl < 2.5(7 away from the mean. This arbitrary cutoff is meant as a compromise 
which shows the effect of nonzero width of the luminosity functions, without extrapolating too 
far into the tails in which there is as yet no data. To give a feel for the sensitivity of our results 
to the assumed luminosity function width we repeated our analysis for luminosity functions with 
larger and narrower \M — Mmeanl ranges, (but with fixed observed intrinsic a). We find that 
the total supernova counts vary less than 0.88% when \M — Mmeanl ranges from 2a to 3cr; this 
insensitivity reflects the fact that the bulk of supernova counts are from events near the means of 
the distribution. One the other hand, we found that the maximum observed supernova redshift (and 
thus the depth to which one can probe the CSNR) is very sensitive to the choice of \M — Mmeanl- 
For example, the LSST maximum supernova redshift in 1 year is ^^s^ — 0.89 for our fiducial 
choice of |M — Mmeanl = 2.5cr, as seen in Table[2j On the other hand for |M — Mmeanl = 2cr and 3a, 
we find 2:sn,max = 0.73 and 1.06, respectively. Here rare, intrinsically bright events determine the 
redshift reach, and the deeper the luminosity function reaches into the bright-end tail, the larger 
the resulting 2:sn,max- Thus we would expect the intri nsically brigh test events, of Type Ibc and 
Type Iln, to give the greatest redshift reach. Indeed, Cooke (j2008l ) has recently illustrated how 



Type Iln events can be mapped out to z > 2 by ground-based 8 meter-class telescopes. 

Of course, all of our forecasts assume that the luminosity functions of each supernova type, 
and the relative frequencies among the supernova types, all remain unchanged at earlier epochs. 
However, it is entirely plausible and even likely that these properties could evolve, e.g., with metal- 
licity and/or environment. These effects are likely to be crucial in determining the true redshift 
reach of future sky surveys, and thus predictions such as ours will improve only as real supernova 
data becomes available with good statistics at ever-increasing redshifts, and one can directly con- 
str ain and/or measure evo lutionary effects. Moreover, the relatively small sample sizes available 
to iRichardson et al.l (l2002l'l could well lead to underestimates of the true range of luminosities of 



each type. For example, (jOezari et al.ll2008l ) very recently report of an unusually bright Type II-L 
event, SN 2008es, with peak magnitude Mb — —22.2, far outside of the absolute magnitude range 
we have adopted for this subtype. 
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Indeed, the enormous statistics gathered by future surveys wiU allow for cross-checks and 
empirical determination of other evolutionary and systematic effects. A major such effect is dust 
obscuration, to which we now turn. 



5. Dust Obscuration: Disentangling the Degeneracies and Probing High-Redshift 

Star-Forming Environments 

The loss of some supernovae due to dust obscuration must be understood accurately and 
quantitatively in order to take full advantage of the large data samples of supernovae which are 
detected. As noted above in ^3.21 currently we have very limited knowledge of supenova extinction 
and particularly its evolution, and most of what is reliably known is based on empirical studies 
of supernova counts. Precisely for this reason, future surveys offer an opportunity to address this 
problem in great detail by leveraging the enormous numbers of supernovae of all types, seen over 
a wide range of redshifts and in a wide range of environments. Here we sketch a procedure for 
recovering this information. 

Future surveys will produce well-populated distributions of supernovae; these encode informa- 
tion about extinction and reddening due to dust. Specifically, in redshift bin Az around z one can 
measure, often with very high statistical accuracy, the apparent magnitude distribution for each 
subtype of core-collapse events. These distributions can be made for all bands, but as we have 
shown, detections and/or light curve information will be most numerous in r and g bands; we will 
focus on these for the purposes of discussion. Within a redshift bin, the distance modulus /i is 
fixed, and the light curve and associated JC-correction should reflect intrinsic variations within the 
core-collapse subtype. 

Thus, for a given core-collapse subtype and redshift z, one can construct histograms of r and 
g peak magnitudes. From redshift and supernova type, one can compute the distance modulus and 
K-correction, and use these to infer, for each event, the dust-obscured peak magnitude M^ust = 
mobs — fJ-iz) — K{z) = Mpeak + ^ where Mpeak is the intrinsic peak magnitude for the event, and A 
is the extinction for this event in its host galaxy. By comparing two passbands we can also evaluate 
colors, for example g — t = -M^^peak ~ Mr peak "I" Ei^g — r), where E{^g — r) is the reddening. In 
general, within a redshift bin we expect the A and E to vary on an event-by-event basis, reflecting 
the properties of dust along the particular line of sight through the particular host galaxy. 



Invaluable insight into these issues comes from the iHatano et al.l (jl998l ) analysis of extinction 
in observation of local supernovae. These authors argue that the data are consistent with very 
strong dependence of extinction with the inclination of the host g a.laxy; this alone guar antees that 
A must vary strongly from event to event even within subtypes. iHatano et al.l (|l998l ) also argue 
that the variation of dust column with galac tic radius al so sugges ts that extinct i on is responsible 
for the paucity of supernovae at small radii. (IShawlll979l ). Finally, IHatano et al.l (|l998l ) also point 
out that core-collapse events are more extincted than Type la events because the la's have a higher 
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scale height and thus are more hkely to occur in less extincted regions. 

On an event-by-event basis, intrinsic hght curve and color evolution are degenerate with dust 
evolution. However, the large sample size may allow for a physically motivated empirical approach 
to lifting this degeneracy. If on theoretical grounds we can assume that at least one core-collapse 
subtype has negligible intrinsic evolution in its lightcurve, then for that subtype M and K are 
effectively known and moreover are constant across events in a particular redshift bin. In this 
case, the apparent magnitude and color distributions can be directly translated into distributions 
of extinction and reddening. By comparing these distributions (or e.g., their means and variances) 
across different redshifts, one directly probes dust evolution. 

Moreover, if one can use one core-collapse subtype as an approximate "standard distribution" 
from which to extract dust properties, one might press further by assuming that other core-collapse 
events will be born in similar environments and thus encounter similar extinction and reddening. 
One can thus use the empirically determined dust evolution to statistically infer the degree of 
intrinsic lightcurve variation in the other core-collapse subtypes. 

If subtype can be firmly established, comparison of magnitude distributions of different core- 
collapse subtypes allows for a purely empirical approach, namely, one can compare the evolution 
of the magnitudes distributions of different subtypes. One could first provisionally treat each 
supernova subtype as a "standard distribution" with no intrinsic evolution; then for each subtype 
one would infer dust extinction and reddening at each redshift. It is reasonable to expect that 
the different subtypes sample the same dust properties, as long as the host environments are 
not systematically diffe rent for the d i ffereii t subtypes (all of which arise in massive-star-forming 



environments). Indeed, iNugent et al.l (|2006l ) have performed such an analysis to use V — I colors 



of Type II-P events to infer reddening for events out to z ~ 0.3. 

A comparison of the dust extinction inferred the different subtypes amounts to a test for 
intrinsic variation. With information from multiple subtypes, it may be possible to isolate dust 
effects common to all, and intrinsic variation peculiar to each subtype. For example, if one subtype 
distribution evolves significantly more than another (e.g., one subtype variance grows more than 
another) then the difference in variance must be intrinsic, and that the lesser variance is an upper 
limit to the variance due to dust effects. 

The ability to empirically measure extinction depends on the intrinsic width of the Ajz) dis - 



tribution, and on surveys' ability to probe this distribution. At low redshift, iHatano et al.l (119981 ) 
find a wide (> V^^^) range of extinctions, much of which they attribute to inclination which will 
remain an issue at higher redshift. On the other hand, as a given survey pushes to higher redshift, 
progressively less of the distribution is observable. For the case of LSST, we see in Fig. [10] that 
with nrv^^ = 23™**^ the least obscured events are visible out to redshift z ~ 1, while those which 
have suffered A^- = 1™^§ of extinction would correspond to the m\f^ = 22™^^ curves, which reach 
to about z ~ 0.5. Thus over this shallower redshift range, extinction can be probed in detail, but 
with a narrowing observable range at higher redshift. 
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If future surveys can empirically determine effects of dust evolution, this would not only remove 
a major "nuisance parameter" for supernova and cosmology science, but also gain information 
of intrinsic interest. Namely, we will learn about the cosmic distribution and evolution of host 
environments of supernovae and thus of star formation. 



6. Discussion 



The large amount of core-collapse supernovae observed by synoptic surveys will yield a wealth 
of data and enormous science returns. Here we sketch some of these. 



6.1. Survey Impact on the Cosmic Supernova and Star Formation Histories 



As we have indicated in the previous section, synoptic surveys will determine cosmic core- 
collapse supernova rate with high precision out to high redshifts. Moreover, with the large number 
of supernova counts, and with light curves and host environments known, the total cosmic core- 
collapse redshift history can be subdivided according to environment and/or supernova type. For 
example, with photometric data alone one can determine to high accuracy correlations between 
supernova rate and host galaxy luminosity and Hubble type. One can compare supernova rates 
in field galaxies versus those in galaxy groups and clusters. Using galaxy morphology one can 
investigate correlations between supernovae and galaxy mergers. With the addition of spectroscopic 
information one can also search for correlations with host galaxy metallicity. 

The CSNR is also tightly related to the cosmic star-formation rate. Therefore with the high 
precision CSNR, and assuming an unchanging initial mass function, one can make a similarly precise 
measure of the cosmic star-formation rate. On the other hand, one can test for environmental 
and/or redshift variations in the initial mass function, by comparing the supernova rates based on 
direct survey counts with the star- formation rates determined via UV and other proxies. 

In addition to core-collapse explosions, synoptic surveys will of course by design also discover 
a similarly huge number of Type la supernovae. Thus the Type la supe rnova rates can be com- 
pared to those for core collapse everits. As has been widely noted (e.g., iGal-Yam k. Maoa |200J; 
Watanabe et al.l Il999l : lOda k. Totanil l2005l : IScannapieco &: BildstenI l2005l . and references therein) 
this will yield information about the distribution of time delays between the core-collapse and ther- 
monuclear events. Moreover, one can explore differences in the environmental correlations for the 
two supernova types (and subtypes). 
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6.2. Survey Supernovae as Distance Indicators: the Expanding Photosphere 

Method 



Type la supernovae have become the premier tool for distance determinations at cosmological 
scales, thanks to their regular light curves, high peak brightnesses, and relatively less dusty environ- 
ments. Nevertheless, given the importance of the cosmic distance scale, and the ongoing need for 
systematic crosschecks and calibration, it is worthwhile to consider other methods. Core-collapse 
events offer just such a method, via the expanding photosphere/expanding atmosphere method. 



This method was originally c onceived by iBaadd (|l926l ) and IWesselinkI (j 19461 ) for study of 
Cephieds; iKirshner &: KwanI (jl974l ) applied the Baade-Wesselink method to supernovae. The key 
to the technique is to exploit the simple kinematics of a newborn supernova remnant: the freely- 
expanding photosphere grows in size as R = vt. Thus, for purely blackbody emission, the luminosity 
grows with size (i.e., time) as L = AnR'^aT'^. With good sampling to measurements time t since 
explosion, and spectroscopic inference of v and T, one can recover the luminosity. In principle, 
therefore, one can use the explosion as a standard candle. 

In practice, this method has been slow to mature. Until recently, th e agreement with inde - 
pendent distance measures has been only good to within a factor ~ 2 (e.g., IVinko &: TakatsI 120071 ). 
The complex (out of local thermodynamic equilibrium) spectra of supernovae has proved difficult to 
adequately model. However, recently important advances have been made in the radia tion transfer 



modeling of young supernova remnants and its fitting to spectra of local supernovae (iBaron et al 



2004 : iDessart &: Hillied l2005l . 120081 ). Because of this, the expanding photosphere (or more prop- 
erly, expanding atmosphere) method now appears to be reaching consistency with other distance 
measures; this method now shows agreement appro a-ching the ^ 10% leve l. Similar pecision now 



seems possiblue using a separate, empirical method (jHamuv &: Pintdl2002l ) which exploits the ob- 
served correlation between luminosity and expansion velocity of Type II-P events. This opens up 



core-collapse supernovae as alternative di stance indicators. Indeed, several group (jNugent et al 



20061 : iPoznanski et al.ll2007l : IOlivaresll2008l ) have already applied this method to various collections 
of Type II-P observations, yielding tight Hubble diagrams out to z ~ 0.3. 

To use this method as it is currently envisioned, follow-up spectroscopy is mandatory for each 
event (see ^6.5p . with photometric surveys identifying the candidates. Obviously, for the largest 
surveys, in practice only a tiny fraction of core-collapse events could be studied in a (separate) 
spectroscopic campaign, particularly given that the most common core-collapse types are intrinsi- 
cally dimmer than Type la events and thus require longer exposures to obtain spectra. Followup 
requirements thus are the limiting factor for the use of core-collapse events as distance indicators. 



The situ ation for Type 



Recent work (Poznanski et al 



Wangj|2007l : iKuznetsova et al 



a supernovae is better-studied and also potentially more hopeful . 



2007 



Kim 



fc Miauell 



20071 : iKunz et al.ll2006l : iBlondin fc Tonrvll2007l : 



20081 : Sako et al. 2008 ) suggests that photometric redshifts of Type la 



events near maximum light could be obtained with sufficient precision (give a low-redshift training 
set) to provide useful dark energy constraints without spectroscopy. Whether photometric-based 
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distances can be derived for core-collapse events with sufficient accuracy remains to be seen. It 
nevertheless seems to us a worthy object of further study. In t his context it is worth noting that PES 
plans to do followup spectroscopy on ~ 25% of Type la events (jThe Dark Energv Survev Collaboration 



20051 ). We suggest that at least some modest fraction of this follow-up time be dedicated to core- 



collapse monitoring. 



6.3. Other Science with Cosmic Supernovae 



The physics, astrophysics, and cosmology of cosmic core-collapse supernovae is a fertile topic; 
with our detectability study in hand, a wide variety of problems present themselves. Here we sketch 
these out; we intend to return to these in future publications. 

The huge harvest of core-collapse events will open new windows onto other aspects of supernova 
physics. For example, the physics of black hole formation in supern o vae, a nd the neutron-star/black 
hole divide, remain important open questions. iBalberg &: Shapiro! (120011 ) have estimated the rates 
of events with observable signatures of black hole formation; LSST should provide a fertile testing 
ground for these predictions. 

The elaboration of the cosmic history and specific sites of high-redshift supernovae will also 
offer unique new information about supernova "ecology" ~ i.e., feedback and cycling of energy, mass, 
and metals into the surrounding environment. For example, large surveys will offer the opportunity 
to study supernova rates as a function of host galaxy and galaxy clustering, shedding new light 
onto large-scale star formation and its connection with galaxy evolution. Moreover, DES and other 
surveys will discover an enormous number of rich galaxy clusters; the occurrence of both Type la 
and core-collapse events in clu sters will offer importan t new insight into the origin of the very high 
metallicity of intracluster gas (jMaoz &: Gal-Yamll2004l : iMaoz et al.ll2005l ). 



Core-collapse supernovae also are the sources, directly or indirectly, of high-energy radiation 
of various kinds. For example, supernovae act as accelerators of cosmic rays. These in turn 



interact with interstellar matter to produce high-energy 7-rays. iPavlidou &: Fields! (12002! ) used 
then-available estimates of the cosmic star-formation rate to show that this 7-ray signal has a 
characteristic feature, and makes a significant part of the extragalactic 7-ray background around ~ 1 
GeV. With the successful launch of the high-energy 7-ray observatory GLAST, this component of 
the 7-ray background may for the first time be clearly identified. Regardlessly, a sharper knowledge 
of the cosmic supernova rate (and thus cosmic-ray injection rate) will work in concert with GLAST 
observations to probe the history of cosmic rays throughout the universe. 
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6.4. Comparison with Type la Survey Requirements 



The characteristic of Type la supernovae are in general very similar to the core-collapse su- 
pernovae. Hence the observational requirements for their identification in synoptic surveys are also 
very similar. The rest-frame, full-width at half-maximum timescale for Type la supernovae is ~ 20 
days. Therefore surveys will need a scan cadence of a few days in order to get a well-sampled light 
curve. For ex ample, the Pan-STARRS strate gy for Type la discovery is to sample the light curve 



every 4 days (|Tonrv et al. 



2003 



Tonryl (2003) • As we have discussed, this sampling frequency is 



also suitable for the core-collapse supernovae which the time scale of the light curve also last a few 
weeks. 



6.5. Redshifts and Typing from Photometry and FoUowup Spectroscopy 

Survey supernovae become scientifically useful only when one can establish, at the very least, 
their redshift and whether they are core collapse or Type la. Since followup spectroscopy will 
not be possible for the large numbers of future events, photometric redshifts will be needed. For 
events in which a host galaxy is clearly visible, one can use photometric redshifts of the hosts. 
Here one is helped by the ability of surveys to stack all of the many (non-supernova) exposures to 
obtain a much deeper image than those with the supernovae. Once the host re dshift is know, the 



supernova type must be det ermined. Baysean analysis techniques and software (jPahlen &: Goobai 



2002 



Poznanski et al.l 120071 ) have been developed to distinguish both Type la and core-collapse 
events. These authors find that type discrimination depends crucially on the accuracy with which 
the redshift is known. For spectroscopic redshifts, their methods is extremely accurate, and for 
photometric redshifts the method is still quite good, though in this case misclassifications can reach 
15 — 25% depending on Uphoto-z- 

Followup spectroscopy on a subset of events will be essential to calibrate the accuracy of the 
photometric typing (and host redshifts). In particular, spectroscopy will be invaluable in identifying 
and quantifying catastrophic failures in the typing algorithms; on the basis of these it may be 
possible to refine the routines. As noted in the previous section, followup is also required for any 
events one hopes to use in distance determinations. 

For events without clear host galaxies and without followup, one must resort to photometric 
redshifts and typing based on the supernova light curve itself, in whatever bands a re available. It is 



not cl ear that this can be done with any reliability on an event-by-event basis As iPoznanski et al 



(|2007l ) emphasize, one might make statistical statements about the types and redshifts of the entire 
class of "hostless" events. Here spectroscopic followup will be essential, not only for determining 
the supernova redshift but also the nature of the underlying host. 
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7. Conclusions and Recommendations for Synopic Surveys 

The next ten years will witness a revolution in our observational knowledge of core-collapse 
supernovae. Synoptic sky surveys will reap an enormous harvest of these events, with tens of 
thousands discovered in the near future, culminating with of order 100,000 seen annually by LSST. 
These data will reveal the supernova distribution in space and time over much of cosmic history. 
The needed observations are naturally a part of the scanning nature of these surveys, and require 
only that core-collapse events be included in the data analysis pipeline. 

The potential science impact of this unprecedented supernova sample is enormous. We have 
discussed ways in which the photometric supernova data alone will contribute in significant and 
unique ways to cosmology and astroparticle physics, as well as to studies of core-collapse and 
supernova evolution themselves. We illustrate one such application by demonstrating how to recover 
the cosmic supernova rate from the redshift distribution of supernova counts in synoptic surveys. 
The large datasets ensure that the statistical error will be very low, and the first large surveys will 
rapidly determine the CSNR to precision exceeding that of current data based on observation of 
massive-star proxies. 

With the addition of spectroscopic followup observations, the survey-identified core-collapse 
supernovae can be used as distance indicators. Thanks to recent advances in the phenomenology 
of supernova spectra and the modelling of their expanding atmospheres, the early light curves 
provide standardizable candles. This expanding photosphere/atmosphere method could provide a 
cross-check for the cosmic distance scale as inferred from Type la supernovae. 

To summarize our recommendations for synoptic surveys, in order to capitalize on this poten- 
tial: 

1. Include core-collapse supernovae (all Type II as well as Types lb and Ic) in the data analysis 
pipeline. 

2. Include a scanning mode in which the depth mf^ is as large as possible, in order to maximize 
the supenova redshift range. Surveys which modes which probe down to mf^^^ = 26™^^ could 
discover many supernovae (both core-collapse and Type la) approaching redshift z ^ 2. 

3. Adopt scanning cadence of revisits every ~ 4 days, in order to capture core-collapse events 
at peak brightness, and to obtain a well-sampled lightcurve. This timescale also appropriate 
for Type la events. 

4. Allocate some followup spectroscopy to core-collapse events. This will calibrate photometric 
Type la/core-collapse typing and typing among core-collapse subtypes, and will be particu- 
larly crucial for probing the nature of events in which a host galaxy is not seen. 

We close by re-emphasizing that these recommendations require only modest efforts in anal- 
ysis, little to no modification of the strategies already in place for Type la searches, and some 
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commitment of followup spectroscopy. Thus a small extra investment of resources will reap hand- 
some scientific rewards as we open our eyes to the incessant rise and fall of these beacons marking 
massive star death throughout the cosmos. 



It is a pleasure to thank Joe Mohr, Vasiliki Pavlidou, Tijana Prodanovic, Jon Thaler, Yun 
Wang, David Weinberg, and Michael Wood-Vasey for valuable discussions. We are particularly 
indebtted to Josh Frieman, Avishay Gal- Yam, John Beacom, and David Branch for detailed con- 
structive feedback on an earlier draft, which has greatly improved this paper. 



Note: In the final stages of writing this paper we became aware of the work of I Young et al 



(j2008l ). These authors discuss core collapse rates in sky surveys, with a focus on events in low- 
metallicity environments. Where it is possible to compare and when we adopt their parameters 
(particularly mj^j^), our analyses seem to be in broad agreement. 



A. Appendix: The Supernova/Star- Formation Connection 



The star- formation rate and supernova rate for any astrophysical site are intimately related. 
Moreover, in many applications such as ours the timescales of interest are much longer than the 
~ few Myr supernova progenitor lifetimes. In this case, the star-formation rate and supernova 
rates are proportional. This is expressed above in eq. ([2]). The constant of proportionality can be 
obtained from the initial mass function ^(m). In stellar mass range [m,m + dm) the number of 
new stars is proportional to ^(m) dm, while the mass of new stars is m^(m) dm. Thus the number 
of supernovae per unit new star mass-i.e., the conversion between star- formation and supernova 
rates-is 

7^gN _ /sN ^("^) _ JsN ^("^) dm /sj^ m i{m) dm _ Xsn 



J m ^(m) dm /sn i{m) dm J m (,{m) dm (m) 



SN 



(Al) 



where = /gj^ m ^ dm / J mS^ dm is the mass fraction of new stars that will go into supernovae, 
and {m)sN = /g^ rn£_{m) dm/ Jgj^Cim) dm is the mean supernova progenitor mass. For illustration, 
consider a Salpeter IMF ^(m) oc m~^'^^ over mass range (O.S MfD, IOOMq) and S,{m) oc m~^'^ for the 
low mass range (O.IMq, O.SMq) (IBaldrv Glazebrookll2003l . their "Salpeter A" mass function). We 
also take supernova progenitors to lie in the mass range mgN G {8Mq,50Mq). This gives Xsn = 
0.15,(m)sN ~ I5.95M0, and thus a star- formation/supernova conversion factor A'sn/("^)sn = 
0.00914Mq^ The uncertainty here is significant, probably about a factor of 2. 

The cosmic star-formation rate can be estimated from a number of observables tied to massive 
(i.e., short-lived thus "ii istantaneous" ) star - formation. Proxies often adopted are the UV and/or 
Ha luminosity densities (iMadau et al.lll996l ). Of these, UV light has a more direct co nnection with 



mass ive stars, but is also affected more by the dust extinction than the Ha light (jStrigari et al 
2OO5I ). Most cosmic star- formation studies find a sharp increase in the rate up to z = 1, but there 
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re mains a larg e unce rtainty of the star formation rate at higher redshift. In this paper, we adopt 
the lCole et ahl (j200ll ) fitting formula for the cosmic star-formation rate 

a + bz 



■cole 



/iMoyr-^Mpc" 



(A2) 



" 1 + {z/cY 

where (a, 6, c, d) = (0.017, 0.13, 3.3, 5. 3), which are one o f the best-fitted parameters according 
to current observing data found by iHopkins &: BeacomI (120061 ) . Using the conversion factor of 
Hopkins &: BeacomI (120061 ) and our adopted Hubble constant, this gives a local rate of 7^sn(0) = 



1.1 X 10 SNIIMpc yr . The benchmark CC SNe rate rises to a peak around z = 2.5 and then 
slowly declines at high redshift. 

To illustrate the impact of different star formation, we also did all calculatio ns with an alter- 



nativ e CSNR. Here we normalize to the current observed, c ounting-based CSNR (IBotticella et al 



20081 ) , and take the shape from the a fitting function of the iHippelein et al.l ([2003) star formation 
rate. We also lower the rate by 30% bec ause we want the alternativ e CSNR to be as much different 
as the benchmark CSNR as possible and lHopkins &: BeacomI ([200a) suggested that the uncertainty 
of the normalization of the star formation rate is about 30%. The differences between this and our 
fiducial rate gives a sense of the current rough but not perfect agreement between the CSNR as 
inferred indirectly from progenitor light (sometimes reprocessed) and directly from counting. 



B. Appendix: Supernova Predictions for Upcoming Synoptic Surveys 



Table 3: Sensitivity to survey depth in r band, for fixed sky coverage. 







DES 






Pan-STARRS 




LSST 


[mag] 


Total SNe 


Redshifts 


Total SNe 


Redshifts 


Total SNe 


Redshifts 


21 


4.32 X 10^ 


0.081 < 


z 


< 0.18 


3.24 X 10^ 


0.007 <z < 0.36 


2.16 X 10^ 


0.008 < z < 0.36 


22 


1.86 X 10^ 


0.066 < 


z 


< 0.38 


1.40 X 10^ 


0.007 <z < 0.56 


9.32 X 10^ 


0.008 <z < 0.56 


23 


6.86 X 10^ 


0.064 < 


z 


< 0.66 


5.14 X 10^ 


0.007 <z < 0.89 


3.43 X 10^ 


0.008 < z < 0.89 


24 


2.19 X 10^ 


0.063 < 


z 


< 1.10 


1.64 X 10^ 


0.007 <z < 1.34 


1.10 X 10^ 


0.008 <z < 1.33 


25 


6.32 X 10^ 


0.063 < 


z 


< 1.62 


4.74 X 10^ 


0.007 <z < 1.88 


3.16 X 10^ 


0.008 <z < 1.88 


26 


1.55 X 10^ 


0.063 < 


z 


< 2.17 


1.16 X lO'^ 


0.007 <z < 2.48 


7.73 X 10^ 


0.008 <z < 2.47 


27 


3.17 X 10^ 


0.063 < 


z 


< 2.68 


2.38 X lO'^ 


0.007 <z < 3.08 


1.58 X lO'^ 


0.008 < z < 3.08 



Figures ITVIIOI show the supernova forecasts for different surveys. For each we show the expected 
annual supernova harvest AA^obs(-2) as a function of redshift, with the sky coverage held fixed to 
the values in Table [TJ The count distribution across redshift bins are directly proportional to the 
differential supernova rate distribution rsN,obs(-z)) via eq. (fTS[) . Thus the shapes the curves follow 
those of rsN,obs(-z) as seen in Figure Ulan explained in the accompanying discussion. 

Figures ITlllOlalso show the survey scan time required for these data to constrain the cosmic star- 
formation rate in each redshift bin to within a statistical precision cr(7^sN)/^SN = 10%. We have 
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seen (eg . [201) that the precision at each bin is inverse with the counts, cr{TZsT<s)/TZsN = ANohs{z)- 
Thus these panels show trends in which monitoring time decreases with the counts per bin. This 
mirrors the behavior shown in Figure [6] and explained in the surrounding discussion. 



Table [3] shows the effect of survey limiting magnitude on redshift range and total supernova 
harvest in the r-band. We see that each unit increase Amf^ = 1^^^ in survey depth yields a large 
enhancement, in the total supernovae seen. The supernova numbers increase by a factor 4.3 when 
going from mf^ = 21™^^ to 22™^^ to factor of 2.0 when going from 26™'^^ to 27™^^. Of course, there 
is a tradeoff in the needed exposure. For the faintest objects at the highest redshifts, background 
noise dominates, and monitoring time grows by a factor (lO^/S)^ 6.3 per magnitude. Thus, 
including a narrower but deeper survey mode will likely yield fewer supernovae, but if judiciously 
implemented, this tradeoff may be worth the additional redshift coverage. 
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Fig. 7. — Sloan Digital Sky Survey (SDSS-II). All the results plotted here are using the benchmark 
CSNR. Top panel: The number of supernovae observed in one scan season of 3 months per year, in 
redshift bins of with Az = 0.1. Results are shown for a fixed scan sky coverage AOgcan = 300 deg^, 
and the survey depth as labeled. Bottem panel: The monitoring time At needed in order to 
determine the cosmic supernova rate to a 10% precision. 
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The i^T- Correct ion 



The X-correction accounts for redshifting of the supernova spectrum across the passbands in 
the observer frame. In the context of the Southern inTermed iate Redshift ESO Supernova Search, 
the elegant and instructive analysis of lBotticella et al.l (j2008l ) determine ii'-corrections for the ~ 90 
supernova confirmed and candidate events in their survey. They found that the corrections depend 
strongly on redshift, light curve phase, and on waveband. In particular, the shifts from observed V 
and R bands to rest-frame -B-band both typically have K < 0, i.e, a negative correction, particularly 
at early times most relevant here; this reflects the blue colors of the early phases. The corrections are 
at early times (within the first three weeks) usually a shift |i^| ^ ]^mag^ with the largest corrections 
K ~ —2^^^. As seen in eq. ([7]), a negative correction would reduce the apparent magnitude and 
thus improve the observability of the supernovae. 



In this paper we compute X-corrections following the formalism of iKim et al.l (|l996l ) in a 
single band, and add the color correction term (which described as iJxb in Section 13. 4p to make the 
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Fig. 8. — As in Figure [3 for the Dark Energy Survey (DES). Here results are show for different 
mjj^, but for a fixed scan sky coverage AOscan = 40 deg^. 
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corresponding correction when transfer from one band to another. In band x, the if-correction is 

fF{X)S^{X)dX 



= 2.5 log(l + z) + 2.5 log 



J F{\/{l + z))S^(X)d\ 



(CI) 



where F{\) is the unobscured, rest-frame spectral distribution of the supernova, and S'x(A) is 
the sensitivity of filter x. We included five types of core collapse supernovae: lb, Ic, II-L, II- 
P, Iln. For each spectral ty pe, we adopt a rest-frame spectrum F{\) following the prescription 
of iDahlen &: FranssonI (|l999l ). who adopt blackbody spectra (sometimes slightly modified) with 
different temperatures and time evolution. We picked the temperatures which last for about a 
weak around the peak luminosity and treat them as a constant. Since the surveys we are interested 
in will have cadence less than a week, this should be a reasonable simplification. For Type lb and 
Ic, we choose a 15,000 K blackbody with cutoff at A < 4000A because of UV blanketing. We choose 
11,000 K for Type II-L, 10, 000 K for Type II-P, and 14,000 K for Type Iln. 

Our resulting X-correction appear in Fig. [TTl plotted for the g and r bands which we will see 
are the optimal for supernova discovery. The huge turnoff of the ii'-correction in Type I-bc is due 
to the short-wavelength cutoff in its spectrum. A large value is set so that there will be no Type 
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Fig. 9.— As in Figure El for Pan-STARRS. We hold fixed A^scan = 30000 deg^. 
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I-bc supernovae observed beyond the cutoff point. Regardless of the turnoff of Type I-bc, for low 
redshifts the ^'-c orrection is in th e rang e of +2™*^^ to — but typically is negative, which is 

in consistent with iBotticella et al.l (|2008l ). 
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Fig. 10. — As in Figure [H for Large Synoptic Survey Telescope (LSST). We hold fixed AQg 
20000 deg^ 
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0.0 0.5 1.0 1.5 2.0 

redshift z 

Fig. 11. — ET-correction of the four types of supernovae(I-bc,II-L,II-P,II-N) in g and r bands (green 
lines for (7-band and red lines for r-band). The sharp upturn in the Ibc correction reflects a cutoff 
in the supernova spectrum due to UV line blanketing; see discussion in text. 
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